Single observable concurrence measurement without simultaneous copies 
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We present a protocol that allows us to obtain the concurrence of any two qubit pure state by 
performing a minimal and optimal tomography of one of the subsystems through measuring a single 
observable of an ancillary four dimensional qudit. An implementation for a system of trapped ions 
is also proposed, which can be achieved with present day experimental techniques. 
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Even if entanglement was spotted as a key feature of 
quantum mechanics since the early days of the theory Q , 
it was only with the advent of quantum information sci- 
ence that a great deal of attention was drawn upon the 
problems of characterizing, properly quantifying, and ul- 
timately measuring entanglement |2|]. 

Until recently, measurements of entanglement had only 
been achieved indirectly by performing measurements on 
several non-commuting observables of the system, and 
then adequately combining the results Q . A direct mea- 
surement of concurrence (previously shown to be a proper 
entanglement measure Q), however, was reported in y|, 
in which two copies of the state were used simultane- 
ously, following the idea in Q . Although this experiment 
constitutes a landmark on the path towards fully under- 
standing quantum entanglement, a simpler measurement, 
which does not involve simultaneous copies, is desirable. 

In the simplest case where we deal with a pure state 
of two qubits, one way to address the problem is to take 
advantage of a well-known relation between the bipartite 
concurrence and the reduced density matrix of either sub- 
system p|: 



C 2 =4det 



Pq, 



(1) 



where p q is the reduced density matrix of one of the 
qubits and C is the bipartite concurrence. Hence, the 
concurrence of the system can be obtained by perform- 
ing the tomography of only one of the qubits. 

In ref. [8( , Rehacek et al. presented a protocol for opti- 
mal minimal qubit tomography, in which all information 
pertaining to the state of one qubit is obtained by mea- 
suring the population of the states of two ancillary qubits, 
which are previously entangled with the target qubit by 
nonlocal operations. In a loose sense, information of the 
measured qubit is "written" in the ancillas: the three val- 
ues (fr x ), (fry) and (fr z ) necessary for qubit tomography 
are encoded into the four probabilities Pjk (j, k = 0, 1) 
of the ancillary system to be found in each state \jk), of 
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which only three are independent because of the unity 
sum requirement. An implementation where the qubit 
was encoded in linear photon polarizations and different 
paths in an optical interferometer played the role of an- 
cillary qubits was also proposed. Later on, this proposal 
was experimentally achieved 9] . 

The aim of this article is twofold: first, we want to 
point out the fact that concurrence for a two qubit pure 
state can be obtained through measurement of the prob- 
ability distribution of the spectrum of a single observ- 
able, without the need of simultaneous copies of the state. 
This is achieved by performing a minimal and optimal to- 
mography of one of the qubits via a single measurement 
on an ancillary four dimensional system. The tomogra- 
phy is minimal in the sense that no redundant informa- 
tion is obtained from the measurements (as opposed to 
the standard procedure) and optimal in that it achieves 
maximum accuracy in determining an unknown state • 
Even if the procedure developed in || could in principle 
be used for such an end, we will use in order to illustrate 
our point a tomographic protocol of our own, in which the 
fact that the probability distribution of the spectrum of 
a single observable is being measured appears naturally. 
By introducing a four level qudit as our ancillary sys- 
tem, we are able to "write" the three desired values in 
the populations of the four levels of the ancilla Pq , Pg> , 
Pe and Pgy, of which only three are independent. By 
choosing these states to have different energies, we can 
pick the observable to be the energy of the ancilla. 

Second, we propose an implementation of the proto- 
col for a system of trapped ions, which is achievable 
with present day experimental techniques. Our proto- 
col proves simpler to implement for this kind of systems 
than that in [8|, for it uses one less ancilla, which means 
that one less ion is involved. 

In what follows, we will denote by \x) the two qubit 
pure state on which tomography of one qubit is to be 
performed, and \G), \G'), \E) and \E') the four distinct 
states of the ancilla. We will make use of two kinds of 
operations: 

i) Rotations between ancillary states. These are de- 
noted by R^ K (9) = exp (— if <r^), where fr£ K is one of 
the Pauli operators (a € {x, y, z}) defined on the sub- 
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space spanned by the arbitrary states | J) and \K): 
\K)(J\ + \J)(K\ 



JK 



a 



a J y K = -i{\K){J\-\J){K\) 
ai K 



(2) 



\K){K\-\J){J\. 



ii) Controlled operations applied on the target qubit 
and controlled by the ancilla. These are denoted by 
C A U, where A <E {G, G', E, E'} denotes the control an- 
cillary state whose occupation implies action of opera- 
tor U on the selected qubit. We can have, for instance: 

C E U^(\G) + \E))\ X ) = ^(\G) + \E)U)\ X ), where 

we explicitly put operator U, which acts only on the tar- 
get qubit, to the right of the ancilla ket. Only three 
instances of U will be actually realized: a x , a y and — <r 2 , 
which are unitary operations on the qubit, whose basis 
states are denoted by \g q ) and |e 9 ). 

Protocol. Our protocol starts with the system in the 
state |G)|x)- To this initial state, we apply three succes- 
sive rotations, R GE {6i), R GG '(0 2 ) and R G ' E '{6 3 ) with 
6>i, 02 and #3 such that 



R GG '(0 2 )\G) 



r g ' e '(0 3 )\g') 



1 
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1 
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V5\G) + \E) 
V3\G) + V2\G') 
(\G')-\E')). 



(3) 



We thus obtain the state 



The next step of the protocol requires us to perform the 
controlled operations C (& y ), C E (a x ) and C E (—& z ), 
ending up with the following state: 

^( |G) + ^| GVs + ^ + -L im )lx>. 

(5) 

Finally we apply the following local tt/2 rotations 
around the y axis on the ancilla: R GE (§), R G ' E ' (§), 

R ° °' (f ) and R v E ' (l) to obtain: 

(\G)Q G - \G')Qg> - \E)Q E + \E')Q E ) |x), (6) 
where we have: 
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(7) 



We then readily calculate the probabilities for the an- 
cilla to be found in each state: 



Pg 



i (l + ^=((**) - (* v ) - <*«» 



(8) 



By adding and subtracting these probabilities, we can 
obtain any mean value (d' a ), which means we have suc- 
cessfully performed the tomography of the qubit just by 
measuring the populations of the energy eigenstates of 
the ancilla. 

We would like to stress that the four values of the prob- 
abilities in Eq.|Hlare the expectation values of the opera- 
tors Q G , Q%>, Q% and Q E i, which constitute a minimal 
and optimal POVM 0] for the one qubit tomography we 
are considering. 

Having performed the tomography of one of the qubits, 
it is straightforward to find the value of the concurrence 
of the bipartite system using relation JIJ. in terms of the 
occupation probabilities, it is given by: 



C 2 = A(l-Z{P 2 G + P*,+P 2 E + P 2 E ,)) 



(9) 



We thus managed to obtain the value of the concur- 
rence by measuring the probability distribution of the 
spectrum of a single observable of a four dimensional an- 
cillary system, with no need of simultaneous copies of the 
state. We stress the fact that even if only the expectation 
value of an observable is needed instead of the complete 
probability distribution in schemes using simultaneous 
copies, in practice this distribution must be determined 
anyway in order to compute the expectation value [5|- 

Trapped Ions Implementation. Consider now a system 
of ions inside a linear Paul trap. To a good approxima- 
tion, the effect of the trap in the motion of the ions can 
be described by a harmonic oscillator. The qubits are en- 
coded in ground and excited electronic states of each ion, 
while one of the ions plays the role of the ancillary system. 
When an ion is illuminated by laser light quasi-resonant 
with one of its electronic transitions, the collective mo- 
tional degrees of freedom can be coupled to the electronic 
ones via photon-momentum exchange. The laser excita- 
tion can be done in several different ways, giving rise to a 
large number of possible interaction Hamiltonians. Here 
we will be interested in a situation where the motional 
sidebands are well resolved and the so-called Lamb-Dicke 
limit applies. Moreover, we will only consider the exci- 
tation of one collective motional degree of freedom, the 
center of mass (CM) motional mode in the longitudinal 
trap direction. 

In order to perform the operations requested by the 
protocol, we have to consider the laser excitation of any 
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given electronic transition of an ion in three different 
ways. One consists in illuminating the ion with laser 
light resonant with the transition, often called carrier ex- 
citation. The second way is to excite the transition with 
light resonant with the first lower motional sideband (red 
sideband); and the last one uses light resonant with the 
first higher sideband (blue sideband). 

Under the conditions stated above, the interaction 
Hamiltonians corresponding to each one of these situ- 
ations are given, in the interaction picture, by [Tlj : 



4P 



H\n\e i,p 5+ +h.c. 
r?fi|Q|e i<?s (7+a + h.c. 



^|ft|e i0 a+ a 1 +h.c, 



(10) 

(11) 
(12) 



respectively. Here the operator <r+ is the electronic rais- 
ing operator, and a and are the annihilation and cre- 
ation operators of the CM vibrational mode, respectively. 
il = |f2|e l< ^ is the laser Rabi frequency and r\ is the Lamb- 
Dicke parameter, which in the Lamb-Dicke limit satisfies 
rj < 1. 

The rotations appearing in the protocol are performed 
via carrier excitation of the ions, according to the time 
evolution operator: 



-j-iicT 



Uc(0,4>) = e 



r(cosO) sin(<j>) CT» e ) 



(13) 



where 6 = |£1|t and a 9 ^ are the electronic Pauli operators 
(see Eq. (J3J)), acting on the subspace spanned by \g) and 
|e), these in turn representing generic ground (|G), \G'), 
\g q )) and excited (\E), \E'), \e q )) states, respectively. In 
particular, the rotations R 9 x e (9) = U c (0, 0) and Rf{6) = 

Uc(0, — 7r/2) are obtained by adjusting the laser phase. 

Controlled operations are achieved through excitation 
of the CM vibrational mode. Both Jaynes-Cummings 
(red sideband, Eq. (|ll|l 'l and anti Jaynes-Cummings 
(blue sideband, Eq. H12|0 interactions are used in the 
protocol. 

For concreteness, we present an implementation using 
40 Ca + ions, which have been used in several experiments 
in Innsbruck [12|. States \g q ) and \e q ) of each qubit can 
be encoded in sublevels m — —1/2 of the 4S' 1 / 2 state and 
m = —1/2 of the 3-D 5 / 2 state, respectively. 

We further use another 40 Ca + ion for our four level 
ancilla. States \G) and \G') can be associated to the 
m = —1/2 and the m = 1/2 sublevels of the 4S , 1 / 2 
state, respectively, while states \E) and \E') can be as- 
sociated to the m — 3/2 and m = —3/2 sublevels of the 
metastable 3D 3 / 2 state (see fig. 

The first part of the protocol consists in preparing 
the ancilla in the superposition given by Eq. |0J. Since 
4Si/2 <-* 3Z?3/2 is a quadrupole transition, it is possible to 
select sublevel transitions with Am = ±2 by controlling 
the angle between the incident laser beam and the direc- 
tion of a weak applied magnetic field, as well as the laser 
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FIG. 1: Encoding of the ancillary ion and relevant transitions 



polarization. These transitions allow us to perform rota- 
tions between specific ground and excited levels without 
disturbing the rest. Rotations R y among the two ground 
state levels \G) and \G') are performed through a carrier 
Raman excitation, which can be achieved through simul- 
taneous off-resonant excitation of the dipole transition 
4jSj/2 4P 3 / 2 by two laser beams of polarizations tt and 
(T + focused on the ancillary ion. 

We start with the system in state |G)|x)|0)„, where 
|0)„ is the ground state of the CM vibrational mode. In 
order to perform the first rotation R y 7E (9i) we apply a 
carrier excitation of the 4 ( S , 1 / 2 3Z? 3 / 2 quadrupole tran- 
sition with Am = 2. This rotation is followed by a carrier 
Raman excitation among levels \G) and |G") as described 
above, implementing R^ G (6* 2 ). We finally perform the 
R G e (#3) rotation applying a carrier excitation of the 
4Si/a 



3D 3 / 2 transition with Am — —2. After that 
series of laser pulses, the electronic state of the ions is 
given by Eq. (0J , with the vibrational CM mode still un- 
occupied. 

Next in the protocol are the controlled operations. 
First, we want to apply a y on \x) controlled by state 
\G'). We apply on the ancilla a 7r-pulse resonant with 
the first blue sideband of the 4/S]/2 <-» 3Z? 3 / 2 transition 
with Am = —2, and with phase = 0. (This choice 
of phase will be maintained for all red and blue detuned 
pulses). This anti Jaynes-Cummings interaction takes 
state |G')|x)|0)„ to |-E")|x)|l)i>, leaving the other states 
unchanged. With this operation we transfer the informa- 
tion about the occupation of state \G') to the motional 
state |1)„. We now apply a v on the qubit, controlled by 
the motional state |l)t,. We achieve this by performing a 
carrier rotation R 9 x qeq (— 7r/2) on the qubit, followed by a 
2-7r-pulse resonant with the first red sideband of the tran- 
sition between level \g q ) and an auxiliary level \e' q ), which 
takes |ff 9 )|x)|l)f back to itself through |e^)|x)|0)„, while 
gaining a minus sign. We then apply a carrier rotation 
R a x " eq {-K/2) on the qubit to obtain state -\E')d y \x)\l)v 
Finally, another 7r-pulsc on the ancilla resonant with the 
first blue sideband of the 45"!^ «-> 3£> 3 / 2 transition will 
bring the state -\E')cr y \x)\l) v to \G')a y \x)\0) v . Notice 
that all other states remain unaffected by these transfor- 
mations. 

An analogous procedure is followed in order to act on 
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the qubit with a x , controlled by state \E): we first ap- 
ply on the ancilla a 7r-pulse resonant with the first red 
sideband of the 45 , 1 / 2 «-> 3Z?3/2 transition with Am = 
2, which transforms state |S)|x)|0)„ into — |G)|x)|l)„. 
Then we apply on the target ion a Ry qeq (n/2) rotation, 
followed by a 27r-pulse resonant with the first red side- 
band of the \g q ) <-> \e' q ) transition, and a R^ e "{— tt/2) 
rotation. Finally, we apply on the ancilla another red 
sideband 7r-pulse, identical to the first one. 

The —a z operation, controlled by level \E') can be 
achieved by applying the following sequence of pulses: 
a 7r-pulse with Am = —2 resonant with the first red 
sideband of the 4S' 1 /2 *-* 3-D 3 / 2 transition on the ancilla, 
a 27r-pulse resonant with the first red sideband of the 
\g q ) <-> \e' q ) transition on the target ion, and another n- 
pulse on the ancilla, identical to the first. After these 
three controlled operations are performed, the system is 
left in the state given by equation J5J) , with no excitations 
in the vibrational mode. 

We complete the protocol by performing the four rota- 
tions: RGB (|), rG'e' r gg' (|) and ^ee' ^ by 

applying successive carrier laser pulses. We obtain thus 
the state given by Eq. @ • 

At this point, we need to measure the populations 
of the ancillary electronic states. This is accom- 
plished, via electronic shelving technique, by exciting 
the iSi/2 *-> 4P| j -x transition and monitoring the fluo- 
rescence light 13] . In our case, a preliminary series of 
laser pulses is needed in order to prepare the ancillary 
ion for measurement. 

First, a carrier Raman 7r-pulse excites the 3D 3 / 2 <-> 
3-D5/2 transition via the 4P 3 / 2 level, using two ir- 
polarized laser beams. This brings the population of 
state \E) to the Zeeman sublevel 3D 5 / 2 (m = 3/2) and 
the population of \E') to 3D 5 / 2 (m = —3/2). Then, a 
7r-pulse with Am = 2 resonant with the 45i/2 «-> 3D 3 /2 
transition is applied, bringing the population of state \G) 
to 3D 3 /2 (m = 3/2), followed by a 7r-pulse with Am = 
resonant with the 4S l 1 / 2 «-> 3D 5 / 2 transition, bringing the 
population of state |G') to (m = 1/2). Finally, a 

7r-pulse with Am = 2 resonant with the 4 1 S , 1 / 2 3-D 3 / 2 
transition is applied, bringing the population of state |G) 
back to level 451/2 (m = -1/2). 

Now transition 451/2 ^ 4Pi/ 2 IS excited and fluo- 
rescence light is monitored. If fluorescence is observed, 
this indicates that the ancillary state \G) was occupied, 
and the measurement ends. Otherwise, a 7r-pulse with 
Am = resonant with the 4S 1 / 2 «-* 3D 5 / 2 transition 
is applied, bringing the population of state \G') back 
to level 45*1/2 ( m = 1/2). The fluorescence test is re- 
peated, a positive result indicating that state \G') was 
occupied. Again, if no fluorescence is observed, a 7r-pulse 
with Am = 2 resonant with the 4S , i/ 2 <-> 3-D 5 / 2 transi- 
tion is applied, bringing the population of state \E) to 
45i/ 2 (m = —1/2). The fluorescence test is repeated 
once more, now a positive result indicating occupation of 
the \E) state. No light observed in this last test indicates 



that state \E') was the one occupied. Several iterations 
of this process yield the occupation probabilities of the 
four ancillary levels. This whole procedure is equivalent 
to measuring the probability distribution of the spectrum 
of a single observable: the electronic energy of the ancilla. 

We wish to note that the whole series of pulses would 
be considerably reduced by taking advantage of the Zee- 
man splitting of the levels obtained with a strong ap- 
plied magnetic field, but this would also limit the initial 
Doppler cooling of the ions 0] . 

We would like to stress that our protocol is designed 
for pure states. However, if we consider small deviations 
from a pure state we still may have a good estimate for 
the concurrence. Assuming, for example, a density ma- 
trix of the form p = Xp' + (1 — A)|x)(xl> where p', is a 
separable state and A« 1, then one can show that the 
difference between the values of the square of the true 
concurrence and the value obtained using the above pro- 
tocol is -2A(1 - P ■ P') + C(A 2 ), where P> and P are 
the Bloch vectors associated with the reduced density 
matrices of p' and |x)(x|) respectively. 

Before concluding we would like to briefly compare 
our protocol with measurements involving simultaneous 
copies Clearly, the main difference lies in the fact 

that we do not require simultaneous copies of the state in 
order to measure concurrence. Another advantage of our 
protocol is that, given that the relation stated in Eq. 
holds for the bipartite concurrence of a qubit with an ar- 
bitrary number of other qubits, provided they are all in 
a pure state, our protocol is trivially generalized to this 
case. It is also important to note that we do not have 
an extra cost for not using simultaneous copies: the four 
dimensions of our ancilla match those in the two qubits 
where the copy of the system is provided. 

As compared to standard tomography, our protocol not 
only involves a single observable, but is also minimal and 
optimal, implying in particular that improved accuracy 
is achieved in determining the state of the system. 

In summary, we presented a minimal and optimal to- 
mographic protocol which involves the measurement of 
a single observable. This in turn is used to obtain the 
concurrence of a pure two qubit state. We also propose 
a realistic implementation of the protocol for a system 
of trapped ions, which could in principle be carried out 
with present day experimental techniques. 
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